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Abstract 

A  method  for  constructing  asymptotic  expansions  of  solutions  of 
the  reduced  wave  equation  (V^  +  k^)u  -  0  is  derived.  The  method  yields 
expansions  containing  exponential  decay  factors  and  fractional  powers 
of  k,  and  the  possible  powers  of  k  in  such  expansions  are  determined. 
The  expansions  contain  a  number  of  undetermined  quantities  which  can 
be  adjusted  to  yield  expansions  of  solutions  of  particular  problems. 
The  expansions  yielded  by  this  method  are  mere  general  than  those  given 
by  the  Luneburg-Kline  method,  and  possess  the  behavior  exhibited  by 
solutions  of  various  diffraction  problems.  In  the  Appendix  a  theorem 
on  the  asymptotic  expansions  of  solutions  of  general  linear  equations 
is  derived. 
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1.  Introduction 

In  many  of  the  diffraction  problems  which  arise  in  electromagnetic  theory. 


acoustics  and  other  fields  it  is  necessary  to  determine  the  behavior  as  k  ->oo 
X 


(k  -  -r-  ,  where  \  is  the  wavelength)  of  solutions  of  the  reduced  wave  equation 


(1)  (7^  *  k^)u  -  0. 

This  behavior  is  usually  determined  by  obtaining  exact  solutions  and  expanding 
them  asymptotically  in  k.  If  it  were  possible  to  obtain  the  asymptotic  expan- 
sion of  a  solution  directly  from  (l)  and  the  other  conditions  of  the  problem  it 
might  be  possible  to  obtain  such  expansions  in  problems  for  which  the  exact  solu- 
tions are  not  known.  It  might  also  be  easier  to  proceed  in  this  way  even  in 
problems  for  which  the  exact  solutions  are  known. 

A  direct  method  for  constructing  the  first  term  in  the  expansion  of  a  period- 
ic electromagnetic  field  was  devised  by  R.  K,  Luneburg "- -^ ,  and  this  method  was  ex- 
tended by  M,  Kline  *•  -'  to  yield  further  terms  in  the  expansion  and  also  extended 
to  other  equations  [llj.  This  expansion  is  related  to  the  rays  of  geometrical 
opticsj  the  first  term  in  the  expansion  is  called  the  "  geometric  optics  field"  . 
Derivations  of  some  of  the  same  results  have  been  given  by  F,  G,  Friedlander'-  -' , 
J.  Riblet  (unpublished),  H,  Bremmer'-  ^   and  E,  T,  Copson'-''i,  This  method  is  ex- 
plained and  applied  to  a  number  of  diffraction  problems  by  J,  B,  Keller  and  I. 
Kay^-*;  another  application  was  given  by  C,  Schensted'-  J. 

The  above  method  does  not  yield  all  asymptotic  expansions  of  electromagnetic 
fields.  Other  ^i^ds  of  expansions  which  contain  exponential  decay  factors  and 
fractional  powers  of  k,  neither  of  which  appears  in  the  Lune burg-Kline  expansions, 
have  occurred  in  various  diffraction  problems.  Such  expansions  have  been  found 
from  the  exact  solutions  for  the  electromagnetic  field  due  to  a  dipole  near  a 
sphere  (H.  Breramer  [7^ )  and  for  the  acoustic  field  produced  near  a  rigid  cylinder 
by  a  plane  wave  (F.  G.  Frie  dlander  \9\)»     Similar  terms  occur  in  an  approximate 
solution  obtained  by  W.  Franz  and  K,  Depperman  [8}  for  the  electromagnetic  field 
near  a  perfectly  conducting  cylinder  due  to  an  incident  plane  wave."*' 
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In  this  paper  we  attempt  to  obtain, by  a  direct  method, more  general  asymptotic 
expansions  than  those  given  by  the  Luneburg-KLine  method,but  we  only  consider  solu- 
tions of  equation  (1),     We  seek  expansions  containing  exponential  decay  factors  and 
fractional  powers  of  k,  in  order  to  include  those  expansions  found  in  the  problems 
referred  to  above,     SpecificaJly  we  take  the   solution  of  (l)  to  be  in  the  form 

(2)  u  =  v.exp(ik</- k°^) 

and  assione  that  v  has  the  asymptotic  expansion 

(3)  V-   21      -^ • 

n"o  ,    n 

k 

The  functions  3.    X     and  v     and  the  real  numbers  a  and  X     are  to  be  determined  in 
'  n  n 

such  a  way  that  \  ^^  >  \     and  that  equations  (2)   and  (3)   satisfy  (l)  formally. 
Conditions  on  all  these  quantities  are  obtained, by  means  of  which  expansions  of 
the  form  given  in  (2),    (3)   can  be  constructed.     However  we  do  not  prove  that  the 
expansion  so  constructed  is  the  expansion  of  some  solution  of  (l). 

In  the  Appendix  to  this  paper  we  state  and  prove  a  useful  theorem  about  the 
asymptotic  expansions  of  solutions  of  linear  equations.     By  means  of  this  theorem 
the  constants  X     in  expansions  of  the  tjrpe  given  in  (3)   can  be  determined,   and 
equations  for  the  v     can  be  deduced^     This  theorem  is  used  in  our  analysis  of  (3), 

2,       Derivation  of  the  Asymptotic  Expansion 

As  a  first  step  toward  the  determination  of  the  expansion  given  by  (2)  and  (3)» 
we  will  show  that  the  constant  a  in  (2)  may,  without  loss  of  generality,  be  assumed 

to  be  positive.     To  this  end  we  observe  that  if  a  is  negative  or  zero,  then  the 

^-/ 
-k  A  ^J 

Taylor  expansion  of  e  in  powers  of  k  ^    proceeds  in  non-increasing  powers  of  k, 

Theirefore  if  this  factor  is  expanded,   and  the  resulting  series  multiplied  by  the  ex- 
pansion in  (3),  the  product  so  obtained  will  also  be  of  the  foi^  given  by  (3),     Thus 

-k  X 
u  can  be  rewritten  with  a  modified  v  and  without  the  factor  e    '       •     But  this  expan- 
sion,  as  well  as  that  obtained  when  a  -  1,  will  coincide  with  the  Luneburg-Kline  ex- 
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pansion, as  can  be   seen  by  applying  the  theorem  in  the  Appendix,     From  that  theorem 

-k  7^ 
it  follows  that  when  e      'is  absent  or  when  a  »  1  all  the  X     are  integers,   as  in 

the  Luneburg-Kline  case.     Therefore  only  positive  values  of  a  different  from  unity 

need  be  considered  in  seeking  new  expansions. 

Inserting  (2)  into  (l)  yields 

(U)       k^[l-(V9)^-  2ik^"^vV9   .  VX  +  ik(2Vv  .  V<^-^yV^(;f\^  k^VvX)^-k"[2Vv.7/+v7^/] 

+  y^v  •  0, 

If  a  >  1,  the  highest  power  of  k  in  (U)  is  2a  and  therefore  its  coefficient  may  be 
equated  to  zeroj  this  yields 

(5)  (vX)^-o. 

We  now  assiane  that   ^  is  real,  or  that  there  exists  a  complex  constant  b  such  that 
b  A    is  real,  and  then  (5)   implies  that    %,    ±e  a  constant,     Ifence  the  factor  e"     ^ 
can  be  removed  from  (2)  and  then  (2),   (3)  again  reduce  to  the  Luneburg-Kline  expan- 
sion.    Therefore  the  case  a  >  1  need  not  be  considered  in  seeking  new  expansions, 
and  a  is  restricted  to  the  range  0  <  a  <  1, 

With  a  so  restricted,  the  first  term  in  (U)  contains  the  highest  power  of  k, 
and  upon  equating  this  term  to  zero  we  obtain 

(6)  {^f  -  1  . 

Equation  (6)  is  the  eiconal  equation  of  geometrical  optics  which  has  been  thoroughly 
analyzed,  and  which  can  be  used  to  determine  the  phase  9,  After  the  first  term  in 
(U)  has  been  eliminated  the  second  term  contains  the  highest  power  of  k;  equating 
this  term,  in  turn,  to  zero  we  have 

(7)  79  .  7  X  -  0. 

Equation  (7)   states  that  the  equal- amplitude  surfaces   /   ■  const,  are  orthogonal 
to  the  equal-phasa  surfaces  ^  *  const.     Therefore  the  surfaces    /L  ■  const,  contain 
the  rays,  which  are  straight  lines  orthogonal  to  the  surfaces  {?  »  const,,   and  hence 

these  surfaces  are  ruled  or  developable. 


-  u  - 

We  will  now  show  that  we  need  not  consider  the  possibility  a  >  ^  ,  for  in 

o  o 

that  case  the  term  k     v(7X)     would  have  the  highest  power  of  k  of  the  remaining 
terms.     Equating  this  tenm  to  zero  would  yield  ecfiation  (5),  and  the  considerations 
immediately  following  that  equation  also  enable  us  to  exclude  this  case.     Thus  a 
is  restricted  to  the  range  0  <  a  <  "t  • 

It  is  now  necessary  to  distinguish  between  the   two  cases  o.  ^  -r  and  a  <  «■  • 
In  the  first  case  the  third  and  fourth  terms  in  (U)  contain  the  same  power  of  k, 
while  in  the  second  case  all  terms  in   (U)  have  distinct  powers  of  k.     Using  (6), 

(7)  in  each  of  these  cases,  we  obtain  instead  of  (U) 

(8)  ik[27v  •   y$?  +  v7^53-ivV^/3  -  ]i}^^{2Vv   .   V/-  +  vV^%3  +  ^^v  -  0,  «  "  f 
and 

(9)  ik[|7v  •   V(jS  *  vV^i2]  +  k^%(vX-)^  -  k°'{2Vv   *   vi  +  vV^X3   +  V^v  -  0,  0  <  a  <  |  . 

We  now  insert  the  expansion  (3)   for  v  into   (8)   and  (9)   and  assume  that  it  can 
be  differentiated  tennwise  to  yield  the  asymptotic  expansions  of  the   derivatives  of 
V,     In  the  resulting  expansions  the  highest  powers  of  k  occur  in  the  terms  resulting 
from  insertion  of  the  first  term  in  (3)   into  the  first  terms  in  (8)   and  of  (9). 
Equating  each  of  these  terms  to  zero  yields 


(10)  2Vv     •   yj?  +  V  V^p  -  iv  V^/   -  0 


a"  J 


(11)  2Vv^   •  vg  +  v^V^^I     -  0  0  <  a  <     i  . 

Doth  equations   (10)   and  (ll)   are  first-order  homogeneous  linear  differential  equations 
for  V  ,   and  they  are   in  fact  ordinary  differential  equations  along  the   rays,  i,«.,  the 
orthogonal  trajectories  of  the  surfaces  ^   =  const,,  which  are  straight  lines.     Equation 
(11)  is  the   same  equation  as  that  satisfied  by  the  first  term  in  the  Luneburg-Kline 
expansion,   and  it  has  a  simple   solution  and  geometrical  interpretation. 

In  order  to  determine  the  constants  K^  and  the  equations  satisfied  by  the  v^, 

for  n  >  0,  we  first  assume  that  no  v     for  n  >  0  satisfies  the  same  eqiation  as  does 
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v^.     If  this  assumption  were  not  made,  then  the  expansion  for  v  could  consist  of  a 
combination  of  any  number  of  expansions  of  different  solutions,  each  starting  with 
an  arbitrary  power  of  k.     With  this  assumption,  we  may  apply  the  theorem  given  in 
the  Appendix;  this  yields  the  following  results: 

a)         A.     is  the  (n+l)-th  of  the  numbers  obtained  by  arranging  the  numbers 
m_(i-2a)  +  ra^Cl-a)  +  m.  in  ascending  order,  where  m,,  m^,  m_  are  any  non-negative 
integers.     In  the   case  a  *  -r  ,   \     '  ^  » 

b)lf  a  '  r;  the  coefficients  v     satisfy 
d  n 

(12)  2Vv^  .   7{?  +  v^(vViV^X)   -  -H^v^  •   y/+  v  V^X], 

(13)  27v^   o   V9  +  v^(7ViV^/  )  =  -d^Vl   *   ^^*  \-l^^/^3   *  ^^^-2'  "  -  ^• 
c)  If  0  <  a  <  ^,  the  coefficients  v     satisfy  the  equations 

(lU)  2W     .  V^  +  v  V^i?  =  i(v/)^v^   ,    .-  i(2Vv     ,    ..  7?^  +  V     ,    vV^X)+  iV^v     ,    ., 

n        '^         n     "^         '   ^   '     P-j^(n)  P2(n)  ^2^^^  PoC")' 

n  >  Oo 
The  integers  pi(n),  P2(n)   and  Po(n)  are  defined  by  the  equations 

(15)  >c  -(l-2a)  -   \     ,    ^,       K  -d-a)   m  \     .  \  -1  "  X     ,    .    , 

n  P^(n)'         n  V2^^^  "  V^M 

If  any  of  these  equations  has  no  solution,  the  corresponding  p. (n)  is  not  defined 

and  the  corresponding  term  is  omitted  from  (lU). 

To  exemplify  these  results  for  the  case  0  <  a  <  t:  ,  we  will  consider  the  case 

a  ■  —  which  occurs  in  the  diffraction  problems  referred  to  in  the  Introduction,  In 

this  case  we  have  l-2a  =  -r,   1-n  »  rr   and  therefore  X  ■  -r  .  Then  from  (l5)  we  see 

3       3  n   3 

that  p,(l)  »  0  while  P2(l)  and  p^(l)  are  undefined}  p-|,(2)  -  1,  P2(2)  -  0  whiie 
p-(2)  is  undefined,  and  for  n  >  3,  P^(n)  »  n-i  (i  »  1,2,3).  Thus  from  (lU)  we 
obtain 
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(16)  2Vv^  •  V{2  +  v^y^(?  -  i(V/^)^  Vp, 

(17)  27V2    •  V^  +  VgV^J?  =  i(v/)^v^  -  i(2Vv^   •  vX+  v^V^/), 

(18)  2Vv^  •  V9  +  v^7^;2  -  i(^^)^Vr^^^^V2*^^*V2^^'^^*  ^^\-3»     "  -  ^* 

1  13 

As  another  example,  let  a  ■  y.     Then  l-2a  ■  p-  ^  l-<i  -  j-  aiid  we  find  that 

X     ■  0,  X     -  -r-  for  n  >  0,     Thus  the  K     form  a  sequence  of  multiples  of  — ,  with 
the  term  -r  itself  omitted.     From  (15)  we  find  that  p,(l)  •  0  while  Pp(l)   and  p,(l) 
are  undefined;  P2(2)   =  0  while  p-,(2)   and  p-(2)   are  undefined;  p,(3)  *  i,  p.^(3)  -  0 
while  P2(3)  is  undefined;  and  for  n  >  3,  P- (n)  •  n-i-1  (i»l,2,3).     The  equation 
for  the  V    may  be  immediately  obtained  from  (lU)  with  the  aid   of  these  results, 

3.       Conclusion 

A  procedure  has  been  found  for  the  construction  of  asymptotic  expansions  of 
the  form  given  by  (2),   (3)  which  formally  satisfy  (1).     According  to  this  procedure 
the  phase  function  9  may  be  any  solution  of  the  eiconal  equation  (5).     The  level 
surface  of  the  amplitude  function  ^  are  generated  by  the  rays,  which  are   straight 
lines  orthogonal  to  the  surfaces  9  "  const,;  and  the  values  of    p^    on  these  surfaces 
may  be  assigned  arbitrarily.     The  constant  a  must  lie  in  the  interval  0  <  a  <  ^  • 
If  a  »  ^  then  X     ■  ^  and  the  v     can  be  determined  successively  from  (10),   (12)   and 
(13)  which  are  linear  ordinary  differential  equations  along  the  rays.     If  0  <  a  <  p 
then  K    is  the  (n+l)-th  of  the  numbers  obtained  by  arranging  the  set  of  numbers 
m-.(l-2a)  +  m2(l-a)  +  ra^  in  an  increasing  sequence,  where  m, ,  m^  and  m,  are  any  non- 
negative  integers.     In  this  case  v     can  be  determined  from  (ll)  and  then  the  other 
V     can  be  obtained  successively  from  (lU),  ail  of  which  are  also  linear  ordinary 
differential  equations  along  the  rays. 

The  arbitrary  elements  in  this  construction  are  the  initial  values  of  ^  and 
of  the  V     on  some  s^arface,  the  values  of  jL  on  each  ray,    and  the  value  of  a. 
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These  cfiantities  may  be  adjusted  in  order  that  the  expansion  correspond  to  the 
solution  of  a  particular  boundary  value  problem. 

Appendix 
Consider  a  linear  operator  L  and  a  solution  v  of  the  equation 

(Al)  Lv  -  0. 

Suppose  that  L  possesses  an  asymptotic  expansion,  or  a  convergent  expansion, 
with  respect  to  a  parameter  e  at  e  «  0,     Thus 

00       a. 

(A2)  ^^11     ^       \t  ^i+l  >  "i  ^  0» 

i»o 

Here  the  L.   are  also  linear  operators  and  the  c.   are  non-negative  real  numbers 

with  a     ■  0«     Then  it  is  possible  that  v  also  possesses  an  asymptotic  expansion 

of  the  form 

00       A. 
(A3)  V  ^    E     s    J,j    ,  Xj^,  >  Xj    , 

here  the  \.  are  undetermined  real  numbers  and  the  v.  are  iindetermined  elements  of 
J  J 

the  same  type  as  v(e.g.,  functions,  vectors,  etc.).  If  v  has  such  an  expansion, 
what  values  may  the  K.   assume,  and  what  equations  must  the  v.  satisfy? 

From  the  linearity  and  homogeneity  of  (Al)  it  is  clear  that  we  may  multiply 
V  by  an  arbitrary  power  of  e  without  affecting  the  equation,  so  obviously  K     is 
arbitrary.  To  obtain  conditions  on  the  other  K.   and  on  the  v,  we  insert  (A2)  and  (A3) 
into  (Al),  assuming  that  the  asymptotic  expansion  of  Lv  is  obtained  by  termwise 
application  of  the  L^  to  the  v..  Thus  we  have 

00    a.+K. 


(Aii)  E   ^  ^  ^Vj  '^°* 


J-0 
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The  smallest  eKp)onent  of  e  in  (Ah)  is  a     +  X     "  X  j  equating  the  coefficient 
of  this  power  of  e  to  zero,  we  have 

(A5)  V^  -  0. 

Now  we  equate  the  coefficient  of  e     to  zero,  obtaining 


(A6)  L^v,,  +  Yl    L^v^  -  0,         k  >  0. 


a.+X^=A 

We  now  assume  that  none  of  the  v.    for  k  >  0   satisfies  (A5),   for  then  it  could  be 
the  first  term  in  a  distinct  asymptotic  expansion.     But  then  the  sum  in  (A6)  must 
contain  at  least  one  tern,  so  there  is  at  least  one  pair  a. ,  X.  that  satisfies 
a.    +  X.  =  X.  .     Since  a.    is  non-negative  by  assumption,   X    >  X.  and  therefore 
k  >  j.     Thus  each  X,    is  the  sum  of  a  preceding  X.  and  some  a. •     As  the  first  X 

^  CO 

is  X  this  immediately  implies  that  every  X  is  of  the  form  X  +2^  m.a,  where 

i=o 
the  m^  are  any  non-negative  integers.  Thus  the  X.  have  been  determined  and  the 

equations  (a5),  (A6)  have  been  found  for  the  v..  The  v.  can  be  found  successively 

from  these  equations  starting  with  v  ,  since  the  equation  for  v,  involves  only 

other  v.  with  j  <  k.  Furthermore  all  of  the  equations  (a6)  are  of  the  same  foim, 

and  only  the  inversion  of  L  is  required  for  their  solution. 

These  results  can  be  collected  into  the  following 

Theorem!  Let  L  be  a  linear  operator  with  the  asymptotic  expansion  (A2)  and  v  be 

a  solution  of  (Al)  with  the  asymptotic  expansion  (A3).  If  the  asymptotic  expansion 

of  Lv  is  given  by  termwise  application  of  the  L^^  to  the  v  and  if  L^v^^  f   0  for  k  >  0, 

then  V  satisfies  (A$)  and  v.  for  k  >  0  satisfies  (A6).  The  constant  X  is  arbitrary 

OK  O 

but  K.    is  the  (k+l)-th  number  in  the  increasing  sequence  formed  from  the  set  of 

00 

ni;imbers  X     +  ^    m.a,   where   the  m.    are   any  non-negative  integers. 
i«l 
In  the  application  of  this  theorem  to  the  problem  considered  in  the  body  of  this 

paper,  the  L^    are  differential  operators,  of  which  there  are  three  in  (8)   and  four 
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in  (9)  and  e  -  k"^. 

It  is  intei*esting,   and  often  convenient,  to  note  another  way  of  deriving  the 
asymptotic  expansion  (A3)   for  v.     To  this  end  we  insert  (A2)   into   (Al)   and  multiply 
on  the  left  by  L^   .     We  then  obtain  a  different  asymptotic    "  equation"    for  v, 
naiteiy 

X 


(A7) 


V  '>j   e       V 


00       a 


Here  v^  is  a  solution  of  (A5)  and  X  is  an  arbitrary  real  number.  If  we  now  re- 
place  V  on  the  right  side  of  (A?)  by  e  v  we  obtain  an  approximate  expression  for 
V.  If  we  then  replace  v  on  the  right  side  of  (A7)  by  this  now  approximate  expression, 
and  repeat  this  process  indefinitely  we  obtain 


(A8) 


00 


00 


v-L     -E 

n»o  V   i=l 


1  T-lr  \n   o 
e   L  L.  }  c  V  . 

0   1/       o 


Vte  will  now  prove  that  the  asymptotic  expansion  (A3)  of  v  is  also  obtained  by 
rearranging  the  terms  in  (a8)  in  increasing  powers  of  6«     We  begin  by  equating 
the  expansion  in  (A3)  to  that  obtained  from  (A8)  which  yields 


(A9) 


e^<v 


X  + 
_  o 


2~  m.a.    00  m 


X     + 


y  in. a." 
^—    i  1 


In  (A9)  the  m.   are  any  non-negative  integers  and  the  summation  is  over  all  sets 
of  m.   satisfying  the  condition  indicated. 

From  (A9)  we  see  that  the  K.  determined  from  (A8)  have  the  same  values  as 
were  previously  found.  It  is  also  clear  that  the  v  satisfy  (A6),  since  if  we 
multiply  (A9)  by  L     and  rearrange  the  resulting  equation  we  have 


(AIO) 


I-  V.   -    E 

ft    Ir  t.— . 

S,3 


a  +X.-  K 

8      j         K 


E 
i™iy 


00 


m. 


7  (-'■A' '  'o 

i»l 


m.a." 
V.W  i  i 


i 
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In  (AlO),  m^  denotes  the  same  set  of  non-negative  integers  as  does  m.  in  (A9) 
except  that  one  positive  m,  that  one  with  i  «  s,  is  reduced  by  unity.  From  (A9) 
we  see  that  the  bracketed  expression  is  just  v.,  so  (AlO)  coincides  with  (A6)j 

J 

this  completes  the  proof  of  the  statement. 
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